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ABSTRACT 

The direct study of molecular gas in inner protoplanetary disks is complicated by uncertainties in the 
spatial distribution of the gas, the time- variability of the source, and the comparison of observations 
across a wide range of wavelengths. Some of these challenges can be mitigated with far-ultraviolet 
spectroscopy. Using new observations obtained with the Hubble Space Telescope-Cosmic Origins Spec- 
trograph, we measure column densities and rovibrational temperatures for CO and H 2 observed on 
the line-of-sight through the AA Tauri circumstellar disk. CO A - X absorption bands are observed 
against the far-UV continuum. The CO absorption is characterized by logio(iV( 12 CO)) = 17.5 ± 0.5 
cm~ 2 and T ro t (CO) = 500^ 2 qq K, although this rotational temperature may underestimate the local 
kinetic temperature of the CO-bearing gas. We also detect 13 CO in absorption with an isotopic ratio 
of ~ 20. We do not observe H 2 absorption against the continuum; however, hot H 2 (v > 0) is detected 
in absorption against the Lya emission line. We measure the column densities in eight individual 
rovibrational states, determining a total logio(N(H2)) = 17.9^0 3 cm~ 2 with a thermal temperature of 

T(H 2 ) = 2500^700 The high-temperature of the molecules, the relatively small H 2 column density, 
and the high-inclination of the AA Tauri disk suggest that the absorbing gas resides in an inner disk 
atmosphere. If the H 2 and CO are co-spatial within a molecular layer ~ 0.6 AU thick, this region is 
characterized by (n# 2 ) ~ 10 5 cm~ 3 with an observed (CO/H 2 ) ratio of ~ 0.4. We also find evidence for 
a departure from a purely thermal H 2 distribution, suggesting that excitation by continuum photons 
and H 2 formation may be altering the level populations in the molecular gas. 
Subject headings: protoplanetary disks — stars: individual (AA Tau) 
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1. INTRODUCTION 

Molecular hydrogen (H 2 ) is the prima ry constituent o f 
giant planets, both in the s olar system flStevensonlll982l) 
and in extrasolar systems (jSudarskv et al.l 120031) . H 2 is 
typically assumed to make up the majority of the mass 
in disks where giant planets form. Owing in part to its 
lack of permanent dipole moment, H 2 does not emit or 
absorb strongly in spectral regions easily accessible from 
the ground, and is therefore challenging to observe di- 
rectly. The second most abundant molecule in these 
disks, carbon monoxide (CO), is often observed as a sur - 
rogate for H 2 (|Salvk et alj|2007t iRodriguez et al.ll2010D . 
The use of CO as a tracer for H 2 , and therefore as a 
tracer of the molecular gas mass of the inner disk, relies 
on a CO/H 2 conversion factor. In practice, this con- 
version factor is only well-constrained obser vationally in 
diffuse and translucent interstellar clouds (Burg h et al.l 
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l2007f ). In dense clouds, model calculations of the CO/H 2 
ratio (defined here as A f (CO)/7V(H 2 )) tend towards 1 - 
3 x 10~ 4 (see e.g., Visser et al. 2009 and references 
therein). This is in agreement with an extrapolation of 
the translucent cloud trend, but has only been verified 
on one sightline (|Lacv et al.lll994D . The CO/H 2 ratio in 
low-mass protoplanetary disks, in particular at planet- 
forming radii (a < 10 AU), is almost completely uncon- 
strained by direct observation. 

A better understanding of the relationship between H 2 
and CO in the inner regions of protoplanetary disks can 
give a more complete picture of the composition and 
structure of planet-forming regions. In addition to ques- 
tions of abundance, the distribution of molecular gas and 
its physical characteristics are of interest for models of 
giant planet formation and their migration (see Najita 
et al. 2007 for a review of molecular gas in inner disks). 
The far-ultraviolet (far-UV) bandpass provides access to 
the strongest dipole-allowed bands systems of both H 2 
and CO. Due to large opacities of dust and gas, far-UV 
photons do not probe the disk midplane directly. How- 
ever, far-UV spectroscopy offers the most direct means of 
simultaneously measuring H 2 and CO in the upper layers 
of inner disks. 

Far-UV spectroscopy has proven to be an effective 
technique for probing the environments of young cir- 
cumstellar disks. Firs t with the Interna t ional Ultravi- 
olet Explorer (IUE) (|Brown et all I198H iValenti et al] 
2000), then with the Far- Ultraviolet Spectroscopic Ex- 
plorer (FUSE; Wilkinson et al. 2002; Herczeg et al. 
2005) and the Hubble Space Telescope- Space Telescope 
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Imaging Spectrograph (HST-STIS) and Goddard High- 
Resolution Spectrograph (GHRS; Ardila et al. 2002; 
Herczeg et al. 2002; Walter et al. 2003), and today with 
£f5T-Cosmic Origins Spectrograph (COS; Ingleby et al. 
2011; France et al. 2011b). The majority of these obser- 
vational efforts have focused on the detection of H2 fluo- 
rescence from hot molecular, inner disk (a < 2 AU; Her- 
czeg et al. 2004) gas illuminate d by Lyq (I Ardila et all 
20021) and o ther atomic lines ([Wilkinson et all 120021: 
France et all 120071 ). The order of magnitude gains in 



spectroscopic sensitivity afforded by COS have made it 
possible to expand these studies to the ultraviolet emis- 
sion lines of CO (France et al. 2011a), enabling direct 
comparison of the dominant molecul ar sp ecies in these 
objects . Recently, lYang et all ([201 1[ ) and iFrance et all 
(|2011al ) have presented the first detections of H 2 and CO 
absorption, respectively, in the UV spectra of Classical T 
Tauri Stars (CTTSs). In this paper, we present unique 
HST-COS observations of a sightline through the upper 
atmosphere of the AA Tauri circumstellar disk where we 
detect the absorption lines of both H2 and CO. 

Our paper is laid out as follows: In §2, we describe the 
COS observations and data reduction. In §3, we discuss 
the qualitative features of the AA Tau spectrum and the 
data analysis performed to derive the H 2 and CO col- 
umn densities in the AA Tau disk. We use these results 
to constrain the physical conditions of both molecules in 
§4, including the spatial distribution of the gas and the 
role of UV photons in regulating the molecular level pop- 
ulations. We present a brief summary of the paper in §5. 
The details of our modeling and the non-detection of H 2 
against the UV continuum of AA Tau are presented in 
an Appendix. 

2. HST-COS OBSERVATIONS OF AA TAU 
2.1. AA Tau 

AA Tau is a well-studied K7 classical T Tauri 
star in the Taurus Molecular Cloud with an age 
of 1-3 Myr dBouvier et al l Il99l iWhite fc Ghezl [20011 
iKraus fc Hillenbrand! I2009h . The mass- accretion rate 
onto AA Tau based on J7-band observations has been 
observed to be in the range 3.3 - 7.1 x 10~ 9 M yr _1 
(Valenti et al. 1993; Gullbring et al. 1998; White 
& Ghez 2001). Both high-re solution dust disk imag- 
ing ([Andrews fc W illiams 200% and scattered light mod- 
els of the photopolarimetric variability of the circumstel- 
lar material indicate a high disk inclination (i ps 75°; 
O'Sullivan et al. 2005, see also Bouvier et al. 1999). 

Because of its high inclination, AA Tau has become 
the prototypical example of dipolar magnetospheric ac- 
cretion onto the star (Bouvier et al. 2007; Donati et al. 
2010). Periodically the magnetospheric accretion column 
passes in our line of sight, producing redshifted absorp- 
tion in Balmer lines. This absorption often coincides 
with eclipses in the photometric light curve, which may 
be caused by a warp in the inner disk caused by the fun- 
nel flow. The eclipse periodicity of 8.22 ± 0.03 days is 
likely similar to the stellar period. 

Due to its relatively high extinction (Ay = 0.74; Gull- 
bring et al. 1998), far-UV H 2 emission was no t detected 
towards A A Tau by IUE ([Valenti et al.ll2000D . As pre- 
dicted by these authors, this emission was below the sen- 
sitivity limit of IUE. AA Tau exhibits a far-UV excess 
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Fig. 1.— The COS far-UV spectrum of AA Tauri. The data are 
binned to one spectral resolution element for display and represen- 
tative error bars are shown in red. In the top panel, we show the 
full dynamic range spectrum, with the strongest atomic emission 
lines labeled in purple. All other strong emission features in the 
1160 - 1660 A bandpass are attributable to H2 fluorescence. The 
bottom panel shows the weak emission and absorption features as 
well as the far-UV continuum. The majority of the weak emission 
lines are fluorescently excited CO and H2. The detected CO A - X 
absorption bands are labeled in green. 



characteristic of young disks (|Inglebv et all [2001 ). The 
higher-resolution UV spectra presented in this work re- 
veal a wealth of molecular (H 2 and CO) and atomic emis- 
sion features consistent wit h activ e accretion of a gas- rich 
disk (§3.1). iNaiita et all ([2009D present ground-based 
mid-IR observations of the molecular inner disk (a < 0.5 
AU), traced by CO fundamental (Av — 1) band emission 
and atomic emission ([Ne II]) from somewhat larger radii. 
Mid-IR spectra from the Spitzer Space Telescope show 
that the AA Tau disk has a rich chemistry of more com- 
plex molecules ( H 2 Q and OH) and org a nics such as HCN 
C 2 H 2 and C0 2 (ICarr fc Naiffal[2008ri . IBethell fc Berrinl 
(2009) use these results to demonstrate that water va- 
por formation may be common in CTTS disks and that 
water may be the dominant reservoir for oxygen in these 
systems. As the most abundant and stable molecules in 
the inner disk, H 2 and CO are the building blocks for 
this complex chemistry. 

2.2. Observations 

AA Tau was observed by COS on 2011 January 06 and 
07 (program ID 11616) during four contiguous space- 
craft orbits. Target acquisition was performed in the 
MIRRORB near-UV imaging mode through the primary 
science aperture. The far-UV spectrum was obtained 
using multiple settings of the G130M and G160M grat- 
ings in order to achieve continuous, moderate resolution 
(Aw ps 17 km s" 1 f\ spectra in the 1133 < A < 1795 A 

7 The COS LSF experiences a wavelength depen- 
dent non-Gaussianity due to the introduction of mid- 
frequency wave-front errors produced by the polish- 
ing errors on the HST primary and se condary mirrors; 
http : //www. stsci . edu/hst/cos/documents/isrs/ 
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bandpass. Two central wavelength settings were used 
in the G130M mode (A1291 and A1327) at two focal 
plane offset positions (FP-POS = 1,4). The G160M spec- 
tra were acquired at three central wavelengths (A1577, 
A1600, A1623) at the default focal plane position (FP- 
POS = 3). In addition to providing continuous wave- 
length coverage across the far-UV bandpass, the use 
of multiple grating settings minimizes the impact from 
fixed pattern noise inherent to the microchannel plate 
detector. A complete description of the COS instrument 
an d on-orbit perfo r mance characteristics can be found 
in lOsterman et alJ (| 2 1 If ) . The one-dimensional spec- 
tra produced by the COS calibration pipeline, CALCOS, 
were aligned and co added using the custom software pro- 
cedure described bv lDanforth et all ()2010[ ). 

3. RESULTS AND ANALYSIS 
3.1. AA Tau Emission Spectrum 

AA Tau displays a rich emission line spectrum across 
the far-UV bandpass (see Figure 1). The majority 
of the lines can be attributed to H2 fluorescence, ex- 
cited by shock-driven H I Lya (see e.g., Herczeg et al. 
2002) . The signature of p hoto-excited CO A - X emis- 
sion ((France et al.ll2011af) is also present. A A Tau dis- 
plays a prominent 1550 - 1650 A excess that is attributed 
to a combination of collisionally excited H2 and CO 
A - X (0 - 1) band emission that may be fluorescently 
pumped through a coinci dence with C IV (|Bergin et al.1 
l200l iFrance et alJl20lTb llah . 

Species tracing high-formation temperature regions 
(T > 8 x 10 4 K) are present in the spectrum, in- 
cluding N V A1239, 1243, C IV A1548, 1550, and 
He II A1640 A. These lines are likel y formed in 
the funnel flows (|Muzerolle et al.l 1200 If) and active 
chromosphere/trans ition region (|Johns- Krull et al]l2000l 
iBouvier et al. 2007). Lower excitation species such as 
Si III A1206, C II A1334, 1335, and Si II A1526, 1533 A do 
not appear to be strong in the spectrum of AA Tau. The 
atomic species and CO fluorescence merit further inves- 
tigation, however in this paper we focus on the molecular 
absorption spectrum through the AA Tau disk. 

In the following subsections, we describe the analysis 
of the H2 and CO absorption feature s, recently identi- 
fied for the first tim e in CTTS spectra (jYang et "alj|2011t 
IFrance et al.ll20lTal ) . We observe H2 absorption lines su- 
perimposed upon the shock-created Lya profile, and we 
describe the reconstruction of this profile and column 
density analysis of the absorbing gas in §3.2. We also 
detect CO A - X absorption against the far-UV contin- 
uum from AA Tau, and in §3.3, we describe the column 
density determination for this gas. 

3.2. H2 Absorption Against the Lya Profile in AA Tau 

H 2 absorption has been observed against both the stel- 
lar continuum a nd transition region emission lines in Her- 
big A e/Be stars (|Roberge et al.l l2001; Martin- Za'i di et all 
2008). Molecular absorption in UV spectra of high- 
inclination CTTSs was beyond the detection limit of 
space-borne instrum ents prior to the installation of COS. 
iHerczeg et al.1 ((20041 ) presented a detailed analysis of H 2 
in the disk of TW Hya, including a reconstruction of the 
Lya line profile and determination of the H2 opacity in 
the disk atmosphere. Their work suggested that H2 ab- 



sorption could in principle be detected against the Lya 
emission line in gas-rich disks. No other disks were ob- 
served by STIS at high enough S /N and resolving power 
to detect these absorptions, however, H2 absorption in 
CTTS systems has recently been de tected against th e 
wings of the strong Lya emission line (jYang et al.ll201ll ) . 
H 2 absorption is observed against the red-wing of the 
Lya profile in AA Tau, shown in Figure 2. We do not 
observe H2 absorption against the far-UV continuum of 
AA Tau, and upper limits on the H2 for this absorption 
can be found in the Appendix. In the following sub- 
sections, we describe their measurement and derive the 
column density of hot H2 in the A A Tau disk. 

3.2.1. Lya Profile Reconstruction 

The H2 transmission spectrum is the ratio of an inci- 
dent Lya profile to the data. We created a simple em- 
pirical fit to the Lya emission line + interstellar absorp- 
tion in order to reproduce the red side of the line profile. 
First, we remove the 1214.5 - 1216.8 A flux, which is 
dominated by geocoronal H I emission. Ha and H/3 spec- 
tra of AA Tau d isplay absorption ch aracteristic of a an 
inner disk wind (jBouvier et al.l 12007). and the blue side 
of the Lya line profile is heavily suppressed; its flux is a 
factor of ~ 10 lower than the red wing of the line. This 
is presumably due to either an accelerating low-density 
outflow or an optically thick but patchy outflow at a sin- 
gle velocity. For the present work, we assume that the 
majority of the Lya emission profile is formed near the 
star and create a baseline on which to measure the H2 
absorption spectrum. 

We found that a single Gaussian emission line com- 
bined with a rest-velocity H I absorption component 
could provide satisfactory fits to the observed Lya spec- 
trum, shown in Figure 2a. The Gaussian emission line is 
characterized by a peak amplitude of 2.65 x 10~ 12 erg 
cm" - 2 s- 1 A^ 1 , a FWHM of 620 km s"\ and a continuum 
level of 9.5 x 10 -16 erg cm~ 2 s _1 A -1 . An interstellar 
H I profile characterized by AT (HI) = 4.3 x 10 20 cm~ 2 
and oh 1 = 10 km s _1 provides a good fit to the data. It is 
interesting to note that this column density corresponds 
to an extinction of Ay = 0.28, assuming Ry = 3.1 and 
the interstellar ga s-to-dust relationsh ip between N(H I) 
and E(B - V) (|Bohlin et al.l 119781 ) . That is a factor 
of ~ 3 lower than the Ay typically assumed for the AA 
Tau sightline, suggesting that a large fraction of the ex- 
tinction derived from optical and near-IR studies arises 
in circumstellar material with a gas/dust ratio less than 
that in the general ISM. The parameters for the emission 
line may be an oversimplification of actual Lya profile in 
AA Tau, but the interstellar values are constrained by 
the inner and outer edges of the observed H I profile and 
are therefore fairly robust. There was a minor residual 
shape to the fit/data H2 transmission spectrum. This 
w 10 % slope was fit by a spline function and removed. 

3.2.2. N(H 2 ,v,J) 

The AA Tau H2 absorption spectrum is displayed in 
Figure 2b. The strongest H2 absorption lines are labeled. 
After dividing out the model Lya emission line, there 
is 2.8 A of spectral region with high enough S/N for H 2 
profile fitting. We created a multi-component H2 fitting 
routine to measure the column density in the nine absorp- 
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Fig. 2. — (a — left) The red wing of the stellar + shock Lya profile of AA Tau. The data are displayed in black with representative 
error bars plotted in red. The Lya model, including a fit to the interstellar H I column density, is shown in blue. The dashed green line 
represents the local Lya radiation field incident on the disk, (b — right) The H2 absorption spectrum after dividing out the Lya model. 
The data are again shown in black with representative error bars shown in red. The best fit H2 absorption spectrum for bu2 = 5 km s 
is shown in green. Prominent absorption features are labeled on the figure. Table 1 lists the best fit H2 column densities for all observed 
lines. 



tion lines into the Lyman (B 1 E+ - X 1 S+) band system 
that contribute significantly to the molecular opacity be- 
tween 1217.48 and 1220.30 41 We create the intrinsic 
line profiles as a function of column density in a given 
rovibrational [v,J] level. These functions are added in 
optical depth space, a transmission curve is created, and 
this function is then convolved w ith the most recent COS 
line-spread function (Kriss 2011) for comparison with the 
absorption spectrum. The best-fit multi-line H2 model 
is found using the MPFIT routine for a fixed bji2 = 5 
km s _1 (see the Appendix for details of the H2 model). 
Initial conditions were determined by first manually fit- 
ting the H2 spectrum, and in order to remove bias intro- 
duced by these choices, a grid of initial conditions were 
searched. This method produced a total of > 6000 indi- 
vidual realizations of the spectral fitting, and the best-fit 
column densities (N(H2,v,J)) were taken from the mean 
of these distributions. 

This multi-component fitting method was mostly in- 
sensitive to the exact initial conditions, except for the 
(0 - 2)R(2) and (2 - 2)R(9) levels, whose wavelengths 
differ by only 0.01 A (X rest = 1219.09 and 1219.10 A, 
respectively). The total column density at this wave- 
length was robust, however the relative columns in the 
two lines was not. We weighted the individual columns 
by the product of the oscillator strength and the rela- 
tive populations of the two levels at a fiducial hot H2 
temperature (T(H.2) = 2500 K). The oscillator strengths 
and populations of the two lines are [fn(2) = 25.5 x 10~ 3 , 
P R{2) = 5.76 x 10- 4 ] and [f m = 31.8 x HP 3 , P R{9) = 
6.24 x 10~ 4 ], respectively. N(2,2) contributes 0.425 and 
iV(2,9) contributes 0.575 to the total 1219.10 A column. 
The (2 - 0)P(17) line is located at the inflection point 

8 We note that there is a possible absorption on the blue side of 
Lya from the B - X (1 - 1) P(ll) line at \ obs = 1212.55 A, but 
the absorption depth is at the RMS noise level of the data at these 
wavelengths. 



TABLE 1 

H2 Absorption Lines observed against Lya. 
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- 1)R(14) 


1218.52 


14 Ql+0.18 
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±<±.zz_ 2g 


0.0255 
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1.04 
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1.56 


(2 


" 2)P(8) 


1219.16 
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i<±.io_ 2g 


0.0214 


8lf 


1.46 


(0 


-2)P(1) 


1219.37 


14 91+0.12 
J-t.»J-_ 32 


0.0149 




1.02 


(2 


- 0)P(17) 


1219.48 


15 60+ 50 


0.0040 


5 4 +23 
°^-22 


1.85 


(3 


-2)R(ll) 


1220.11 


16 34+°- 78 


0.0213 


qn+15 
9Z -18 


1.80 



a Transitions arc for the B X Y^ - X X Y^ H2 band system. 
b Column density and equivalent width calculations as- 
sume — 5 km s _1 , see §3.2.2. 

c Column densities, oscillator strengths, equivalent 
widths, and energy levels refer to the rovibrational level 
of the ground electronic state out of which the Lya pho- 
tons are absorbed. 

between the Lya line and the far-UV continuum, compli- 
cating its line-shape. This line was fit manually (iV(0,17) 
= 4 x 10 15 cm -2 ) and this value was fixed during the 
fitting procedure to reduce run time. The best-fit H2 
column densities are presented in Table 1. 

The error bars presented in Table 1 were determined 
from the minimum and maximum deviations from the 
mean column densities of the best-fit distribution. In 
cases where the fitting uncertainty errors were less than 
the \-a error bars on the data (shown in Figure 2b), we 
took the errors to be the range of column densities that 
could be accommodated within the uncertainty on the 
data. This approach gives a conservative estimate of the 
uncertainties on the measured column density. Signif- 
icantly larger 6-values would imply turbulent velocities 
greater than the isothermal sound speed for T < 4000 
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K, however, in order to fully understand the parameter 
space, we also determined the best-fit parameters assum- 
ing bn2 — 10 km s . As expected, increasing b reduces 
the best-fit column densities. We found that the column 
densities were typically 0.1 - 0.3 dex smaller at high b. 
Two exceptions to this trend were the columns of the (2 - 
2)P(8) (which increased by « 0.6 dex) and (3 - 2)R(11) 
lines (which decreased by ss 1.3 dex). The equivalent 
widths were calculated from the best-fit column densi- 
ties, assuming that the observed H2 lines are located on 
the "flat part" of the curve-of-growth, and are displayed 
in Table 1. 



3.2.3. Thermal H2 Population 

The measured column densities, normalized by the de- 
generacy of the state (gj), are plotted as a function of 
excitation temperature (T exc = E(v, J)/kg, where ks 
is the Boltzmann constant) in Figure 3. We label the 
absorbing levels in the ground electronic state as [v,J] 
and the rovibrational levels of the electronic state 

as [v ,J 0- The energy levels for the ind i vidual rovibra- 
tional levels were taken from iDabrow ski (1981). Gas in 
thermal equilibrium is well-described by a single Boltz- 
mann distribution, which would appear as a straight line 
in Figure 3, while deviations from this thermal distri- 
bution indicate that non-thermal processes are partially 
responsible for the level populations in the H2 gas. Non- 
equilibrium photo-excitation and excess energy associ- 
ated with H2 formation on grain surfaces are generally 
assumed to be the dominant n on-thermal H2 excitation 
processes in interstellar clouds ( Spitz er fc Cochranpl973|) 
and protoplanetary disks (|Nomura fc Millar! 120051 ) . and 
it seems plausible that both are at work in the molecular 
disk of AA Tau (see §4.1). 

The H2 population is roughly consistent with a ther- 
mal distribution for T exc < 2 x 10 4 K. Assuming that 
the lowest energy states observed in absorption represent 
a thermal H2 population, these lines are used to deter- 
mine the total column density of this gas. Figure 3 shows 
three thermal distribution curves that correspond to the 
best-fit, upper and lower limits to the thermal H2 (the 
solid pink, dashed blue, and dashed green lines), normal- 
ized to the measured column in the [v,J] = [2,1] state. It 
may be that the [2,1] state itself is not thermally popu- 
lated, but in the absence of lower-excitation constraints 
(most of the lines that would constrain the T exc < 5000 
K portion of the excitation diagram are located at wave- 
lengths shorter than the COS G130M band), we choose 
the lowest observed energy level for the normalization. 
With this caveat about the normalization in mind, we 
find that the thermal population can be characterized by 
T(H 2 ) = 2500±7°o K. These values are consistent with 
the ther mal origin of th e Lyq- pumped H2 emission pro- 
posed by iHerczeg et al.l ()2004l ) . The total column den- 
sity of the thermal H 2 is thus logio^Ha)) = 17.86±g:|? 
cm~ 2 . We will present a discussion of the possible evi- 
dence for non-thermal H2 level populations in §4.1. 



3.3. CO Absorption 
While H2 absorptions from the ground vibrational 
state (v - 0) are located at A < 1110 A, the A - X ab- 
sorption system of CO spans the COS far-UV bandpass, 
extending blueward from the (0 - 0) band at A ~ 1544 A. 

We identify nine CO bands from v = - 8 in the COS 
spectra of AA Tau (see Figure 1). Spectral blending 
with photo-excited H2 as well as emission from Si IV 
A1394 and C IV A1548 limit the number of bands that 
can be used for spectral fittin j 117 ! . and the low flux lev- 
els (Fx ^ 1 x 10~ 15 erg cm~ 2 s _1 A -1 ) also make de- 
tailed modeling challenging. Broad CO A - X absorp- 
tion lines were first ident ified in the UV spectrum of HN 
Tau ()France et al.ll2011al) , and we employ a similar tech- 
nique for fitting the CO absorption spectrum of AA Tau. 
To date, all CTTS systems that show A - X absorption 
have very broad CO profiles due to blending of closely- 
spaced rotational lines (suggesting CO temperatures of 
several hundred K) and relatively high disk inclinations, 
arguing that the molecules are located in a warm inner 
disk. 

Each CO band was extracted from the spectrum and 
a low-order polynomial was fit to the most uncontam- 
inated nearby portions of the spectrum for continuum 
normalization. Individual band velocities relative to rest 
were determined to correct for any relative offsets in- 
troduced by wavelength calibration errors or in the co- 
addition process. Profile fitting was used to determine 
the relevant parameters: N( 12 CO), N( 13 CO), T rot , and 
Doppler b- value. Because the thermal width of CO is w 4 
times smaller than that of H2, we restricted our fitting 
to bco < 2 km s _1 . 

CO absorption band profiles were modeled using ro- 
vibrational line wavelengths and oscillator strengths from 
Eidelsberg (private communication) for 12 CO. For 13 CO, 
wavelengths were determined from the ground state en- 
ergies calculated from the mass-independent Dunham co- 
efficients of George et al. (1994) and the A 1 II energy 
levels of Haridass & Huber (1994), while the oscillator 
strengths used were from Eidelsberg et al. (1999). The 
wavelengths and oscillator strengths for the perturba- 
tions, which show up more prominently at higher column 
densities, were taken from Eidelsberg & Rostas (2003). 
In the fitting process, we assumed the ground states for 
both 12 CO and 13 CO were populated following a single 
rotational temperature, T rot (CO). 

Of all the bands, the (1-0), (2-0), (3-0), (4-0), (7-0) and 
(8-0) were the least contaminated by other spectral fea- 
tures so we compared these bands simultaneously with a 
grid of absorption line models. Best-fit parameters were 
determined using a chi-squared metric, limiting the fit- 
ting to normalized flux values below 1.1; this minimized 
the diluting effect the strong emission features would 
have on the chi-squared value. Our analysis finds that 
the CO absorption is characterized by logio(^V( 12 CO)) 
= 17.5 ± 0.5 cm" 2 , logi (7V( 13 CO)) = 16.2 ± 1.0 cm" 2 , 
b CO = 1.2 ± 0.4 km s" 1 , and T rot (CO) = 500±|gg K. The 
best-fit CO model for several bands is shown overplotted 
in orange on the COS spectrum in Figure 4. 



9 The rovibrational levels of the ground electronic state following 
fluorescence would then be [v ,J ]. 



10 IFrance et al.l (1201 lap suggested that C IV may be a photo- 
excitation source for CO in these systems, and that possibility will 
be explored in a future work. 
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Fig. 3. — H2 excitation diagram based on measured rovibra- 
tional state column densities (Table 1). All absorption lines are 
for transitions to the Lyman band (73 - X). Several high- J levels 
at T e xc > 2 X 10 4 K show evidence for non-thermal excitation. 
Thermal population curves (normalized to the column density in 
the [v = 2, J = 1] level) for rovibrational excitation temperatures 
corresponding to the best-fit, minimum, and maximum values are 
shown in pink, green, and blue respectively. These curves imply to- 
tal thcrmalized column densities of 7V(H2) = 7.2, 30, and 3.5 X 10 17 
cm -2 for the three curves. 



The large errors in our best-fit parameters are indica- 
tive of the difficulties inherent in fitting the CO absorp- 
tion in data such as these. Several systematic effects con- 
tribute to the error. Due to the large number of emis- 
sion features, finding suitable uncontaminated spectral 
regions for continuum normalization was difficult. The 
(2 - 0) band shown in Figure 4 (upper right) provides 
an example of the challenges associated with continuum 
placement in CTTSs. These emission features, coupled 
with the low signal-to-noise, also enhance the error in 
our temperature fit; as rotational excitation temperature 
increases, the profile near the low- lying J levels, where 
the lines may be experiencing more saturation effects, 
changes less than at higher-lying J levels, which lie at 
longer wavelengths. The result is a larger error bar to 
higher temperatures. At the lower (v - 0) bands, the 
13 CO tends to overlap with the 12 CO and as such tends 
to complicate the fit; however, for the (7-0) and (8-0) 
bands, the 13 CO is cleanly separated and detected. 

3.3.1. Isotopic Fractionation 

That 13 CO absorption is required to fit the ob- 
served absorption profiles is not unexpected consid- 
ering the high 12 CO column densit ies observed in 
A A Tau (jSonnentrucker et al.1 12007! ) . The best- 
fit logi (iV( 13 CO)) = 16.2 ± 1.0 cm" 2 , implies a 
12 CO/ 13 CO ratio of ~ 20, though we caution 
that the uncertainties on both quantities are large 
( 12 CO/ 13 CO - 20t^°) . This ratio is smaller than 
the local interstellar 12 CO/ 13 CO of w 70 (Sheffer et al. 
2007; see also the review by Wilson & Rood 1994) and 
the isotopic ratios derived from IR observations of young 
stellar objects ( 12 CO/ 13 CO « 100; Smith et al. 2009), 



although the discrepancies are not statistically signifi- 
cant. We briefly consider the primary processes that are 
expected to influence the isotopic fraction in molecular 
environments. 

Isoto pic charge exchange enhanc es the population of 
13 CO ([van Dishoeck fc Blackl Il988! ) . becoming efficient 
at kinetic temperatures < 35 K. While the CO band 
widths strongly rule out T rot (CO) < 100 K, turbu- 
lent mixing in the disk (e.g., Ilger et al. 2004) could 
bring 13 CO formed closer to the midplane to the up- 
per disk where we observe the absorbing gas. 13 CO 
is also selectively photodissociat ed by UV photons at 
the surface of molecular clouds (|Ballv fc Langerl fl982[ ) 
due to the faster onset of 12 CO self-shielding, and this 
process enhances the relative population of 12 CO in 
these environments. While the far-UV ra diation field 
at the surface of a CTTS disk is intense (jFogel et al.l 
12011b iFrance et al.l 1201 lbl ) , the absorption bands with 
the largest 13 C O dissociation probabilities are located 
at A < 1060 A (|van Dishoeck fc Blackl 1 1 9881 ). where the 
continuum is weak. 

A 12 CO/ 13 CO ratio of ~ 20 is marginally consis- 
te nt with lowest inner disk carbon fractions computed 
bv I Woods fc Willacy! (|2009f ). Their models suggest that 
our observed 12 CO/ 13 CO ratios would be expected closer 
to the disk midplane (for the reasons noted above). 
Therefore, if low fractional ratios are found in the warmer 
upper disk, then this material was likely formed in cooler, 
shielded regions and has subsequently been vertically 
transported through the disk. It i s also intriguing to 
note that the models presented by iVisser et al.1 (2009) 
predict isotopic fractions < 30 for 7V(H2) > 10 21 cm" 2 . 
Assuming a typical interstellar CO/H2 ratio (10 -4 ), 
N(H.2) ^ 10 21 cm~ 2 is approximately the expected H2 
column associated with our measured value of iV(CO) = 
10 17 - 5 cm" 2 (see §4.3). 

4. DISCUSSION 
4.1. N on- Thermal H2 Population 

These data most likely probe a sightline through the 
upper atmosphere of the AA Tauri disk, a region where 
non-equilibri u m proc esses are expected to be important. 
IFrance et al.l ()2011aD suggested that non-thermal popu- 
lations could be relevant for both the CO and H2 popu- 
lations in inner CTTS disks (see also the protoplanetary 
disk H 2 models of Nomura & Millar 2005). Their work 
was based primarily on model-dependent emission line 
analyses and analogy to the better-studied c ase of molec- 
ular ex citation in photodissociation regions. Ardi la et al.l 
(2002), using fluorescently excited H2 emission lines to 
reconstruct the ground state level populations in a sam- 
ple of five CTTSs, found that the H2 was considerably 
out of thermal equilibrium. However, more recent analy- 
ses based on HST-STIS observations of CTTSs suggest 
that these levels are consisten t with a thermal popula- 
tion (jHerczeg et ah 2004, 2006). It should be emphasized 
that all of the studies cited above deal with molecular 
emission lines, and it is not yet clear whether the CO 
and H2 emission are tracing the same parcels of gas as 
the CO and H2 absorption. 

Using a direct measurement of H2 and CO absorption 
lines along the same line-of-sight and in the same ob- 
servation, we have a much more straightforward obser- 
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vational basis for understanding the molecular processes 
operating in protoplanetary disk atmospheres. The devi- 
ations from the thermal distributions in Figure 3 could be 
direct evidence for non-thermal excitation in the molecu- 
lar disk atmosphere. It is premature to draw broad con- 
clusions until this behavior is observed in similar systems 
because we only detect three states with strongly non- 
thermal behavior, and two of these (2 - 0)P(17) and (3 - 
2)R(11) are measured in regions of low S/N. However, 
the excitation diagram presented in Figure 3 is intrigu- 
ing, and we briefly describe possible mechanisms for this 
excitation in the following paragraphs. It is also interest- 
ing to note that we only measure significant deviations 
from a thermal population at excitation temperatures 
T exc > 2 x 10 4 K, approxim ately the same energy cut- 
off where lArdila et all (|2002l ) saw the largest non-thermal 
excursions. 

There are two main mechanisms that can drive the 
H2 level populations out of thermal equilibrium^: ex- 
cess energy associated with H2 formation on dust grains 
( "formation pumping" ) and excitation by UV photons 
at a rate faster than the molecules can decay via rovi- 
brational emission lines or collisions ("multiple pump- 
ing"). Initial evidence for a non-thermal component to 
the CT TS H 2 population was observed bv lHerczeg et al.l 
(2002), who found H 2 emission lines pumped out of states 
near the collisional dissociation threshold of H2 in STIS 
spectra of TW Hya. These states cannot be thermally 
populated, and a local source of highly excited H2 ([v, J] 
= [5,18]) is required if these levels persist. Dust grains 
are widely believed to be the pri mary site of H2 forma- 
tion in the present-day Univer se ([Hollenbach fc Salpeterl 
[19701 ICazaux fe T iclcns 20]S) and are abundant in pro- 
toplanetary disks (see reviews by Natta et al. 2007 and 
Dullcmond & Monnier 2010). Numerous authors have 
presented H2 formation rates in molecular clouds (e.g., 
Jura 1975; Black & van Dishoeck 1987; Draine & Bertoldi 
1996). Protoplanetary disk models that include molecu- 
lar formation (e.g., Nomura et al. 2007) have considered 
the effects of grain growth, which can create significant 
deviati on from the propertie s of typical interstellar dust 
grains (jVasvunin et al.lf2011h . However, these authors do 
not explicitly describe the dis tribution of newly f ormed 
molecules in the inner disk. iDraine fe Bertoldil |l996) 
find typical formation states at high rovibrational levels 
([v, J}} ~ [5,9], with (Ef orm ) « 2.9 cV. This may suggest 
that the high energy excursions from a thermal distribu- 
tion we observe in Figure 3 are the beginning of a higher 
energy population of H2 recently formed in the disk it- 
self. Additional calculations of the spatial distributions 
of newly formed H2 in CTTS disks would be useful. 

A calculation of H2 photo-excitation in the AA Tau 
disk is more straightforward, and here we provide a sim- 
ple example of the importance of this process. The rel- 
ative importance of UV line and continuum emission to 
the photoexcitation of trace molecules (CN, HCN, H2O) 
has been examined in some detail (e.g., Bergin et al. 
2003). We argue that UV photo-excitation of H2 can be 

11 A third, possibly speculative mechanism for producing highly 
non-thermal H2 populations is excess energy associated with the 
dissociation of more complex molecules. Rovibrationally excited 
H2 and CO have recently been observed in the spectra of comets, 
consistent with these molecules being the photodissociation prod- 
ucts of formaldehyde (H 2 CO; Fcldman et al. 2009). 



important even in the absence of a strong Lya radiation 
field or a population of T(R 2 ) > 2000 K molecular gas. 
The far-UV continuum produced in the accretion flow 
and magnetically active stellar atmosphere provides suf- 
ficient pumping flux to photo-excite even a cold molecu- 
lar gas. Far-UV continuum spectra can now be routinely 
measured i n CTTS systems with COS, and following the 
analyses of iFrance et al.l (|2011bl laf). we can fit the 912 - 
2000 A continuum spectrum from AA Tau. We find that 
the integrated far-UV continuum radiation field is G aa 
= 10 63 G , where G D is the local UV interstellar radia- 
tion field (G = 1.6 x 10" 3 erg cm" 2 s -1 ; Habing 1968). 
This corresponds to a far-UV continuum luminosity of 
Lftfy 9 x 10 30 erg s" 1 . 

Consider an individual H 2 absorption line that 
will have a large population fraction even at low- 
temperatures (T(H 2 ) < 500 K), the Lyman band (4 - 
0)R(1) absorption line at 1049.96 A. The UV photon 
absorption rate, C\ (in units of photons s _1 ), can be 
calculated by 

C x = h d\ a H2 (1) 

where I\ is the incident continuum photon flux 
(J1050 ~ 2 x 10 10 photons cm" 2 s _1 A" 1 at 1 AU for 
AA Tau), dX is the equivalent width of the (4 - 0)R(1) 
line, and <jh 2 is the line-center absorption cross-section 
of an individual molecule, given by 

Off 2 = — — , — A,: fi (2) 
m e c b H 2 

(|McCandlissl l2003t iCartwright fc Drapatzl 11970ft where 
A, = 1049.96 A and is the oscillator strength 
(ra 0.0155). a H2 of the (4 - 0)R(1) transition is 
1.2 x 10~ 14 cm -2 , which is typical of strong H 2 absorp- 
tion lines at A < 1100 A for b H 2 — 2 km s~ 4 . dX > 0.1 A 
for all temperatures when Af(H 2 ) > 10 18 cm -2 . 

We can incorporate these values to further simplify 
Eqn 1 for the case of AA Tau: 

Clo50 . (2 , , 10 - 5) (i^) 2 (_«*_,) (3) 

where r is the radial distance (in AU) of the H2 from 
the star. To first order, when C\ is larger than the tran- 
sition probabilities (A) of the rovibrational H2 emission 
lines, one expects multiple UV photo-excitations before 
the gas can relax via line emission. The average tran- 
sition probability for the quadrupolar H 2 emission lines 
in the 1 - 4 /xm bandpass (Afjfc) - 5 x 10~ 7 s _1 . The 
pure rotational mid-IR emission lines are another 1-3 
orders of magnitude slower. Therefore, we see that > 40 
far-UV continuum photons are absorbed for every one IR 
emission, in the absence of collisional de-excitation. 

While it is clear that UV photons are important in 
the level populations of disk H 2 , it is not clear that they 
are the dominant factor in this balance. The critical 
density is the ratio of the radiative lifetime of a given 
state (the Einstein A value, with units of s _1 ) and the 
collision rate for de-excitation out of that state, T (in 
units of cm 3 s" 1 ). The critical densities for most rovi- 
br ational states of H 2 ar e < 10 4 cm -3 at T > 2000 
K (|Mandv fe Martinlll993l ). H 2 + H collision rates (T H2 ) 
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Fig. 4. — Select CO A - X (v - 0) absorption bands in the 
AA Tau spectra. The data are plotted in black and representa- 
tive error bars are displayed in red. These bands were fit simul- 
taneously and the best-fit model (§3.3) is shown overplotted in 
orange. The model parameters are logio(./V( 12 CO)) = 17.5 cm" 2 , 
logio(iV( 13 CO)) = 16.2 cm" 2 , T(CO) = 500 K, and b C o = 1-2 km 
s _1 . The high temperature required to fit the data argues strongly 
for an inner disk origin for the absorbing CO gas. 



for gas in LTE in the range 2000 - 4500 K are on the or- 
der of 10" 12 - 10" 11 cm 3 s" 1 for the (1 - 0)S(j") near-IR 
rovibrat ional transitions that ca n be observed from the 
ground (jMandv fe Martini [19931 . Therefore at an aver- 
age density of (n# 2 ) ~ 10 5 cm" 3 (§4.2), the collisional 
de-excitation transition rate is Tji2{nij a ) ~ 10~ 7 - 10 -6 
s _1 for this gas, somewhat less than our estimates of the 
photoexcitation rate by continuum photons at 1 AU. In 
a high-density region such as the inner disk atmosphere 
of a CTTS, strong UV pumping will compete with colli- 
sional excitation and de-excitation to produce level pop- 
ulations th at appear neither exactly thermal nor purely 
fluorescent (jDraine fc Bertoldilll996l ). 

The observations we present in Figure 3 agree with 
t he qualitative predictions of t he irradiated d i sk mo dels 
of lNomura fc Millarl (|2005[) and lNomura et al.l (|2007f l. In 
these models, the UV radiation field thermalizes the level 
populations up to a certain T exc through strong photo- 
electric heating. Beyond this thermalization energy, the 
level populations are determined by the UV pumping 
rates for the relevant transitions. In our COS obser- 
vations of AA Tau, we observe a higher thermalization 
threshold (T exc > 2 x 10 4 K vs. T exc < 1 x 10 4 K) than 
predicted by Nomura et al. (2007; Figs. 12 and 13); 
however, this can most likely be attributed to the radial 
distances where these distributions are evaluated. The 
model calculations are evaluated at 30 - 50 AU, while the 
H2 absorption we observe in AA Tau may be occurring 
in the inner disk (§4.2). At smaller radial distances, the 
more intense UV radiation field will enhance the photo- 
electric heating rate, thermalizing the H2 population to 
higher T exc . 

4.2. The Spatial Distribution of the Absorbing Gas 



We have measured the radial velocities of the absorp- 
tion and emission components of both molecules. We 
find that the H2 absorption lines are shifted by a uniform 
+29 ± 8 km s , slightly high er than the pho t osphe ric ra- 
dial velocity measurement of IBouvier et al.l ([2007) (+17 
km s _1 ) and consistent with the atmospheric Balmer- 
series absorption velocity of +34 km s _1 from lJovl (|1949f) . 
A sample of 10 H2 emission lines located both short- 
ward and longward of the Lya emission line show an 
H2 emission radial velocity of +18 ± 8 km s _1 , consis- 
tent with both the narrow optical line velocity (+4 - 12 



km s" 1 ) of IBouvier et al.l (|2007l) and the Ha emission 



velocity (+13 km s" 1 ) from lJovl (|1949f ) . Due to the spec- 
tral overlap of the CO lines, the radial velocities of this 
species are not as well constrained. Furthermore, emis- 
sion lines coincident with the CO absorption bandheads 
complicate this sharp edge in the data. The average CO 
absorption radial velocity of the v — 1, 2, 3, 4, and 7 
bands is +20 ± 13 km s" 1 . The CO emission radial 
velocity, determined from the bandhead of the A - X 
(14 - 4) complex is +22 ± 19 km s _1 . Therefore, we 
conclude that there is no evidence for high- velocity infall 
or outflow in the molecular material we observe towards 
AA Tau. 

The thermalized H2 and CO populations have sig- 
nificantly different excitation temperatures; however, if 
these absorption lines are both formed in the disk, then 
they are most likely in close spatial proximity, if not 
actually co-spatial. We can put a rough limit on the 
inner radius of the absorbing gas using the line-widths 
of the fluorescent H2 emission. We fit the 12 strongest 
H2 emission lines in the 1260 - 1300 A bandpass, find- 
ing (FWHM) =53+12 km s -1 , which corresponds to 
ri n 0.16 AU (assuming Keplerian rotation and ignoring 
thermal and turbulent contr ibutions to the line-w idth) 
for a 0.53 M central star (IGullbring et al.l I1998[) an d 
a disk inclination of 75° (| Andrews fc Williams! 120071) . 
Taking the half-width at zero intensity (HWZI) of the 
line moves ri n to ss 0.11 AU. An approximate outer ra- 
dius of the absorbing material could be where the lower 
bound on the CO temperature (300 K) is in equilib- 
rium with the radiation field, r out ~ 0.79 AU, assuming 
T» = 4000 K and R* = 1.9 R (from the A A Tau SED 
fitting of O'Sullivan et al. 2005). While these inner and 
outer radii probably do not exactly represent the regions 
subtended by the absorbing molecules, they provide an 
order-of-magnitude estimate of the molecular region. 

The average values of the CO/H2 ratio and H2 density 
in the region probed by our absorption line spectroscopy 
are then: (CO/H 2 ) = N(CO)/N(R 2 ) ~ 0.4 [I 2 ; J] and 

8.4 x 10 4 cm" 3 . Assuming 



N(R 2 )/(r out - r v 



an H2 ortho-para ratio of 3, the total ortho- + para-H2 
collision rate, summed over all possible low er levels, is 
"Yot = 2.85 x 10~ 10 cm 3 s" 1 for J = 9 (|Yang et all 



120101 ) . The CO critical density for J < 9 at 300 K is 
nn 2 ~ 1.8 x 10 5 cm" 3 . Therefore, it is possible that the 
high-J CO populations that drive the determination of 
Trot (CO) are significantly sub-thermal. While the high- 
J levels are populated by collisions, these levels can ra- 
diatively depopulate fast enough that the populations 
are not representative of the local kinetic temperature. 
To explore this possibility, we computed RADEX mod- 
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els (Ivan der Tak et aHl2Q0l using the on-line interfac43- 
For the values of T(H 2 ), AT(CO), b C o, and (njj 2 ) found 
above, we find that CO levels J > 9 are characterized 
by temperatures less than the 2500 K kinetic tempera- 
ture. For 9 < J < 38, the CO excitation temperature is 
483 ± 129 K, in agreement with the T rot (CO) ~ 500 K 
derived in §3.3. Therefore, the CO and H 2 we observe 
in absorption are consistent with an origin in the same 
2500 K gas. 

The average H2 number density derived above assumes 
crude limits to the spatial distribution of the absorbing 
H2 and CO, it is instructive to also consider the limit- 
ing cases. At the high-density limit, one might imagine 
that the entire column of absorbing H2 is in the form of 
a proto gas-giant planet. Taking (r out - fj n ) = 1 d„ we 
see that ng, ~ 3 x 10 6 cm -3 . At the low-density limit, 
we might imagine that the H2 uniformly spans the en- 
tire planet-forming region of the AA Tau disk, or (r out - 
fin) = 10 AU. In this limit, we find njj 2 ~5x 10 3 cm" 3 
and estimate a total H 2 mas^^l of Mh2 ~ 2.4 x 10 22 g, 
which is roughly 4 x 10~ 6 M®. This calculation confirms 
our earlier speculation that these line-of-sight measure- 
ments do not probe the bulk of the molecular distribution 
available for planet formation, and are only sampling the 
tenuous upper regions of the protoplanetary disk. 

4.2.1. Relation to CO Fundamental Emission 

As noted in §2.1, CO emission from A A Tau has 
been detected in fundamental band emission near 
5 fim (Carr & Najita 2008; Najita et al. 2009; 
and see also Salyk et al. 2011). Our UV 
absorption-derived CO column densities and tempera- 
tures are roughly consistent with those derived from 
the mid-IR emission lin es logi (iV(CO)) = 1 7.7±g"| and 
T(CO) = 900 ± 100 K (|Carr fc Naiitall200l . This sug- 
gests that both obser vations may b e tra cing a similar 
molecular population. INaiita et al.l (|2009D find CO line 
widths of FWHM = 145 km s" 1 , suggesting that the 
emission arises from gas in Keplerian rotation within 0.5 
AU of the central star (and likely with a significant con- 
tribution from gas within 0.1 AU). Due to the pencil- 
beam nature of our absorption line measurements, we 
do not have information regarding the velocity broad- 
ening of the CO. The CO emission line widths are ~ 3 
times those found for the UV emission lines of H2 de- 
scribed in the previous subsection, meaning that the CO 
emission is likely originating interior to the H2 in AA 
Tau. N(CO) in the range 10 17 " 18 cm" 2 should be as- 
sociated with H 2 column densities iV(H 2 ) > 10 21 cm -2 , 
several orders of magnitude larger than we observe. This 
discrepancy raises questions about both the physical ori- 
gin of the molecular gas and its composition, and these 
points are addressed in the following subsection. 

4.3. CO/H2 in the AA Tau Circumstellar Disk 

We find that the average observed CO/H 2 ratio in 
the AA Tau inner disk is ~ O.4+Q3. We imagine two 
scenarios that could account for this result. In the 
first, the absorbing CO resides in a relatively low-density 

12 http: //www. sron. rug.nl/~vdtak/radex/radex.php 

13 Assuming z = 0.1 r, the volume of this disk is then 
[(4tt/30) x (r 3 ^ - r?J]. 



medium (uh2 iS 2 x 10 5 cm" 3 ), where the higher ro- 
tational states that drive the CO temperature fits are 
sub-thermal. In this case, the kinetic temperature of 
the gas is ~ 2500 K, the CO and H2 are cospatial, and 
CO/H2 ~ 0.4. A CO/H 2 ratio in the range 0.1 - 1 is sig- 
nificantly larger than the typical interstellar dense cloud 
value (~ 10" 4 ; Lacy et al. 1994 and see Liszt et al. 2010 
for a review of CO/H2 in various phases of the molec- 
ular ISM). The collapse of a dense cloud precedes the 
formation of the disk and protostar, thus finding CO/H2 
> three orders of magnitude above this value suggests 
significant chemical evolution has taken place in the AA 
Tau disk. CO/H2 valu e s of ~ unity were recently sug- 
gested by iFrance et all (|2011aft in CTTS inner disk at- 
mospheres, using an orthogonal analysis comparing CO 
emission to typical iV(H2) values from the literature. 

In the second scenario, the H2 and CO abundance ra- 
tios are interstellar, but the absorbing H2 is spatially 
separate from the absorbing CO. Any CO in the "hot 
H2 region" falls below our detection limit and any H2 
in the "warm CO region" is too cool to create measur- 
able opacity at Lya. An implication of this reasoning 
is that there is a reservoir of warm H2 associated with 
the CO absorption, with typical dense cloud abundances. 
Taking the best-fit CO column density (logi (iV(CO)) = 
17.5 cm -2 ), and assuming a molecular fraction (/h^) of 
0.67 and CO/H2 = 10" 4 , we expect a neutral hydrogen 
column density of logio(iV(H)) = 21.5 cm" 2 co-spatial 
with the absorbing CO gas. The neutral hydrogen opac- 
ity associated with 10 21 ' 5 cm" 2 would completely extin- 
guish the flux from the star to > ± 5 A from the Lya 
line center, which is inconsistent with the Lya profile 
presented in Figure 2a. While Jh2 = 0.67 is ruled out, 
this scenario is tenable for Jh2 jej 0.9. It seems plausible 
that the molecular fraction is high in a warm (~ 500K) 
medium with N{R 2 ) > 10 21 cm" 2 . iBurgh et all l|2010f ) 
present a direct comparison between the CO/H2 ratio 
and fu2 across a range of diffuse and translucent inter- 
stellar sightlines. Extrapolating their Figure 2, one sees 
that the CO/H2 « 10~ 4 is reached at f H2 ~ 0.9. While 
it would be speculative to draw larger conclusions based 
on analogy to cold interstellar clouds, it seems that a 
spatially stratified medium is consistent with our obser- 
vations. 

Is this vast reservoir of warm H2 observable? At ^(Ha) 
< 500 K, the upper rovibrational states of the molecule 
are not appreciably populated. There is not enough 
opacity in the traditional HST far-UV bandpass (A > 
1150 A) to create a detectable signal. At A < 1120 A the 
H2 opacity should rise rapidly due to large populations in 
the low-lying rovibrational levels of the ground electronic 
state, producing an observable absorption signature that 
could directly distinguish between the two scenarios laid 
out above. Additionally, at A < 1120 A far fewer H2 emis- 
sion l ines pumped by Lya are present (|Wood fc Karovskal 
2004), and only one H 2 line pumped by O VI contributes 
significantljQ The decrease in spectral confusion may 
offset the decrease in source signal and instrumental ef- 
fective area across this bandpass. With the installation 

14 The C-X (1-3) Q(3) A1119.08 line may be strong, however 
the branching ratio of the C — X (1 - 2) Q(3) A1075.03 line is small 
(< 0.004) 
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of COS, HST can n ow observe at wavelengt hs as short 
as the Lyman limit ([McCandliss et all 120 10), and new 
medium resoluti on modes (such as the G130M A1222) are 
in development (|Osterman et alJl2010f ). Heavily damped 
line profiles require neither high S/N nor high resolution 
for a robust column density determination, and we sug- 
gest that S/N as low as 2 per pixel may be sufficient 
when the data is rebinned to one spectral resolution el- 
ement. Therefore, we predict that it may be possible to 
directly search for warm H2 in the AA Tau disk with 
future far-UV observations. 

5. SUMMARY 

We have presented new spectroscopic observations of 
the AA Tau protoplanetary disk system. We directly 
observe both H2 and CO absorption lines along the line- 
of-sight through a CTTS disk for the first time. These 
Hubble Space Telescope observations allow us to constrain 
the molecular column densities, relative abundances, and 
temperatures of the two main molecular constituents of 
the disk. We find that the H2 absorption seen imposed 
upon the stellar Lya profile has both thermal and non- 
thermal components, with the thermal material charac- 
terized by logio(-/V(H 2 )) « 17.9 cm -2 and a tempera- 
ture of T(H 2 ) ~ 2500 K. The CO absorption implies 



logi (AT(CO)) w 17.5 cm" 2 and T rot (CO) - 500 K. 
Taken together, the observed CO/H2 ratio of ~ 0.4 may 
either be a direct measure of local abundances or indi- 
cate spatially stratified hot and warm molecular layers. 
Additional modeling of UV photo-excitation and in situ 
molecular formation will be very useful for understand- 
ing the physical state of the low-density molecular gas in 
these systems. 
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APPENDIX 

ABSORPTION MODEL AND NON-DETECTION OF H 2 AGAINST THE FAR-UV CONTINUUM OF AA TAU 

In this Appendix, we describe the creation of an H2 model and use it to place limits on absorption against the far-UV 
continuum of AA Tau. The physical characteristics of the H2 population in the disk can be constrained by comparing 
model H2 absorption spectra with the data. Assuming that the H2 level populations are controlled by collisional 
processes ()Herczeg et al.ll2~004f) . synthetic absorption line spectra are constructed given values of the total U2 column 
density (A(H 2 )), the rovibrational excitation temperature (T(H 2 )), and Doppler broadening parameter (6^2)- While 
in principle H2 absorption lines can be searched for from the short-wavelength cut-off of our observations (A ~ 1133 A) 
to the end of the significant rovibrational population (w 1400 A for T(H2) < 3000 K), practical considerations make 
the 1225 - 1300 A band preferred for constraining H2 absorption against the continuum of an accreting pre-main 
sequence star. At A < 1200 A, ci rcumstellar and inte rstellar reddening are the most severe, the continuum used as 
a background source is declining (|France et al.l 1201 lbl ) . and the reflectivity of the HST Optical Telescope Assembly 
(OTA) is falling precipitously. All of these factors conspire to create low S/N at these wavelengths. At A > 1300 A, 
the rovibrational levels [v,J] are not significantly populated, thus only weak absorptions are predicted. 

The 1225 - 1300 A band is at the peak of the COS effective area and therefore provides the best combination of 
sensitivity and rovibrational population in which to search for an H2 absorption signal. The primary limiting factors 
for this work are the S/N of the continuum and spectral overlap with th e numerous emission features from hot gas 
lines (e.g., N V A 1 239, 1243) and H 2 photo-excited by Lya and O VI (jWilkinson et al.l l2002t [Herczeg et all 120061: 
iFrance et alJl2011bh. We created a bsorption spectra in the 900 - 1400 A wavelength range using the H 2 ools optical 
depth templates |McCandlissl l2003'l for the identification and analysis of H2 lines in the COS bandpass. 

H2 absorption is not detected against the far-UV continuum of AA Tau. In order to constrain the column density of 
the H2 , we define an upper limit as the column density which produces absorption line depths greater than the 1 a error 
bars on the continuum. In the 1225 - 1300 A wavelength range, iV(H 2 ) < 1 x 10 19 cm" 2 for T(H 2 ) > 2000 K. At T(H 2 ) 
< 2000 K, we do not expect a significant number of detectable absorbers redward of Lya, therefore this analysis is not 
sensitive to a high-column, low-temperature H2 component in the outer disk (see §4.3). At T(H.2) > 2000 K, additional 
absorbers are predicted that are not observed, and the increased population in the high [v,J] states lowers the upper 
limit on iV(H2). The 2000 K limit also implies bji2 > 4 km s^ 1 because this is the thermal width in the absence 
of turbulent velocity broadening. In fact, if turbulence is significant, then the column density limits become smaller 
because increasing bn2 increases the absorption depth at line-center for a given A^H^). This result is due to the fact 
that while these lines do not appear to reach zero ("become black") at their line centers when convolved with the COS 
LSF, the strongest lines are indeed saturated. As bji2 increases, the equivalent width of the saturated line increases, 
producing less line-center transmission for a given column density. We see then that our limit logio(iV(H2)) < 19.0 
cm~ 2 is rather conservative because higher temperatures and turbulent velocities only act to reduce this value. 
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